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Abstract. Resonant interactions between ions and Alfven/ion-cyclotron (A/IC) waves may play an important role in the 
heating and acceleration of the fast solar wind. Although such interactions have been studied extensively for "parallel" waves, 
whose wave vectors k are aligned with the background magnetic field Bo, much less is known about interactions between 
ions and oblique A/IC waves, for which the angle between k and Bo is nonzero. In this paper, we present new numerical 
results on resonant cyclotron interactions between protons and oblique A/IC waves in collisionless low-beta plasmas such as 
the solar corona. We find that if some mechanism generates oblique high-frequency A/IC waves, then these waves initially 
modify the proton distribution function in such a way that it becomes unstable to parallel waves. Parallel waves are then 
amplified to the point that they dominate the wave energy at the large parallel wave numbers at which the waves resonate 
with the particles. Pitch-angle scattering by these waves then causes the plasma to evolve towards a state in which the 
proton distribution is constant along a particular set of nested "scattering surfaces" in velocity space, whose shapes have 
been calculated previously. As the distribution function approaches this state, the imaginary part of the frequency of parallel 
A/IC waves drops continuously towards zero, but oblique waves continue to undergo cyclotron damping while simultaneously 
causing protons to diffuse across these kinetic shells to higher energies. We conclude that oblique A/IC waves can be more 
effective at heating protons than parallel A/IC waves, because for oblique waves the plasma does not relax towards a state in 
which proton damping of oblique A/IC waves ceases. 
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INTRODUCTION 



WAVE-PARTICLE INTERACTIONS 



Resonant interactions with Alfven/ion-cyclotron (A/IC) 
waves are a possible mechanism for ion heating in the 
solar corona, solar flares, and the solar wind. Cyclotron 
heating in low-j3 plasmas primarily increases a particle's 
thermal motions perpendicular to the background mag- 
netic field Bo llU], and may thus be able to explain the 
observed temperature anisotropies of minor ions in the 
solar corona [2] and protons in the fast solar wind JH]. 
(Here, j3 = Snp/B^, where p is the plasma pressure.) In 
solar flares, the magnetic tension in reconnected mag- 
netic field lines leads to large-scale flows that can gener- 
ate waves and turbulence. Wave energy is then transfered 
from large scales to small scales by nonlinear wave-wave 
interactions 0,Ht]. Small-scale A/IC waves may be suf- 
ficiently energetic in solar flares to stochastically accel- 
erate ions to high energies Most previous studies 
of ion heating by A/IC waves have focused on "parallel 
waves," for which the angle 9 between the wave vector k 
and Bo is zero. On the other hand, in the solar corona 
and solar flares, A/IC wave intensities are not restricted 
to = 0. In this paper, we thus focus on resonant in- 
teractions between protons and oblique A/IC waves, for 
which 0^0. 



We consider A/IC waves in a low-j3, proton-electron 
plasma, and assume that the real part of the wave fre- 
quency, co/tr, is given by the cold-plasma A/IC dispersion 
relation M, 
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Q. p is the proton cyclotron fre- 
and v'a = Bo/ ^/A-npo is the Alfven 
speed. Protons strongly interact with such waves only 
when the resonance condition, 
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is satisfied, where V|| (vjj is the component of the particle 
velocity v parallel (perpendicular) to Bq, h\ (k±) is the 
component of k parallel (perpendicular) to Bo, and n is 
any integer JMSl- The strongest interaction occurs for 
n=l H. 

Figure [TJ plots w(k n ) for 0=0 (solid line), as well 
as two dashed lines corresponding to 1 +k\\v\\/Q. p for 



1 

0.8 
0.6 
0.4 
0.2 

0, 



V -- ' — — 






\ x 

\/ 1 +k v /n 

A T il " p 




/ \ where v = - 


-0.05 v A 


/ \ 




/ \ 




/ \ 




/ \ 




/ 1 t k H V H /£2 P 




/ where v = 


-0.2 v A 

A 







FIGURE 1. The intersections between lines satisfy the reso- 
nant condition n=l. 



two different values of the proton parallel velocity V|. 
The intersections of these two lines and the w(k„) curve 
correspond to solutions of equation (01 with n = 1 . When 
| vii | <C Va, equations (fTJ and © imply that JH 
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where £ n || res is the value of k^v^/Q. p that satisfies equa- 
tion (O for 0=0, and v p h(v||) is the parallel phase ve- 
locity co/kii of the resonant waves at 9 = 0. 

Resonant interactions between particles and waves 
cause particles to diffuse in the vp — v± plane. Parti- 
cles interacting with a particular wave with wave vec- 
tor k and frequency CO diffuse within the V|| — v± plane 
along a curve for which the particle energy is conserved 
in a frame moving with velocity (o/k\\ in the direction 
of Bo [8, 9] (wave pitch-angle scattering). If protons in- 
teract with A/IC waves with a broad range of k\\ values 
and 0=0, then particles will diffuse along closed con- 
tours in the vp — vj_ plane. These contours are defined by 
the equation rj = constant, where iflll [l2tl 
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If resonant wave-particle interactions control the evolu- 
tion of the proton distribution function /, then / becomes 
constant on surfaces of constant 77. Once / = /(tj), pro- 
tons stop gaining or losing energy from interacting with 
parallel A/IC waves, and parallel A/IC waves are nei- 
ther damped nor amplified (fl [T3I1 . The shape of the con- 
tour 77 = 0.07 5 v\ is shown in the top panel of Figure [2] 
When protons interact with oblique A/IC waves with a 
single nonzero value of 9, they also undergo wave pitch- 
angle scattering along a set of nested, closed contours 
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FIGURE 2. Upper panel: The proton distribution is taken to be 
constant along the 77 = constant scattering contours for waves with 9 = 
0, and to decrease as one moves to contours that are farther from the 
origin. The solid line is the 77 = 0.075v^ contour. Protons interacting 
with oblique waves will diffuse upward along the dashed line, gaining 
energy, and damping the oblique waves. Bottom panel: The proton 
distribution function is now taken to be constant along the scattering 
contours for oblique waves with some nonzero 0, and to decrease as 
one moves to contours that are farther from the origin. The solid line 
illustrates one such contour. Protons interacting with waves with = 
will diffuse down the density gradient along the short-dashed line, 
losing energy, and amplifying the 8 = waves. The long-dashed line 
is a contour of constant energy in the plasma frame. 



in the V|| — v± plane. However, at a fixed vii <C va, the 
parallel phase velocity G)/k\\ of resonant oblique waves, 
denoted v p h(vii ,9), is greater than the parallel phase ve- 
locity of resonant parallel waves, v pn (v|| ) 111 211 . Thus, at a 
fixed point in the vii — v± plane, the oblique- wave scat- 
tering contour has a larger slope than the parallel-wave 
scattering contour, as illustrated in Figure [2] which is 
adapted from 111 211 . In the top panel of Figure [2 we take 
the proton distribution function to be constant along the 
77 = constant scattering contours of 9 = waves, and 
the solid-line curve corresponds to 77 = 0.075v^. Pro- 
tons at V|| = vrj that are scattered by oblique A/IC waves 
with some nonzero value of 9 will scatter along the 
dashed-line trajectory, which corresponds to constant en- 
ergy as measured in a reference frame moving at velocity 
v ph( v 0i 9) along the magnetic field. [We have artificially 
increased v p h(vrj, 9) relative to v pn (vo,0) in both panels 
of Figure |2]to make the figure easier to read.] If we take 
/ to be a decreasing function of 77, then there will be a net 
diffusive flux of protons upward along this dashed line, 
resulting in an increase in particle energy and damping 
of oblique waves. Elsewhere, we have calculated analyt- 



ically the damping rate of oblique A/IC waves assuming 
that/ = /(Tj)G2. 

In the bottom panel of Figure |2l we take the proton 
distribution to be constant along the scattering contours 
corresponding to waves with a single nonzero value of 0. 
One of these contours is now drawn with a solid line. 
Protons at V|| = vo interacting with A/IC waves with 
0=0 will scatter along the short-dashed line in this 
panel, which locally corresponds to an 77 = constant 
curve. If we take the proton distribution to decrease as 
one moves to closed (solid-line) contours that are farther 
from the origin, then there will be a net diffusive flux 
of protons downward along the short-dashed line in the 
bottom panel of Figure [2] In this case, the protons will 
lose energy, and waves with 0=0 will be amplified. 

Based on these arguments, we make the following 
conjecture. If some mechanism generates high-frequency 
A/IC waves with a range of values, and if the form and 
evolution of the proton distribution function are dom- 
inated by wave-particle interactions, then interactions 
involving waves with nonzero will act to make the 
constant-/ contours steeper in the V|| — v± plane than the 
77 = constant scattering contours of the parallel waves. 
This in turn will lead to the amplification of waves with 
0=0 and cause the angular distribution of the waves 
at large k„ that resonate with the protons to become 
sharply peaked around = 0. Wave-particle interactions 
will then become dominated by waves with 0=0,/ will 
become approximately constant along surfaces of con- 
stant r/, and oblique waves will be damped. In the next 
section, we describe numerical calculations that support 
this conjecture. 

NUMERICAL CALCULATIONS 

In the quasilinear theory of resonant wave-particle inter- 
actions, protons diffuse in velocity space as described by 
the equation 

-zf = hm — V / d^k- ; Gv_l 
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where G= (l—k\\v\\/(Dfr)d/dv± + (k\\v±/(Dfr)d/dv\\, 

Yn,k=E^J n+l (k 1 y_L/^p)+E^J n -i(ki_vx/^ P ) (where 
we have set E kz — 0), V is the volume, J n is the Bessel 
function of order n, Err = E kx ± iEty, and E# (B k ) is 
the Fourier transform of the electric (magnetic) field [8, 
lit]. To integrate equation © numerically, we discretize 
velocity space using cylindrical coordinates with 100 
grid cells spanning the interval < v± < 0.5 va and 100 
grid cells for the interval —0.5 va < V|| < 0. We discretize 
k-space using spherical coordinates with 30 grid cells 
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FIGURE 3. The power spectrum of the smallest-6* waves are 
shown at different times, and are compared to the fixed power 
spectrum of the waves with 6 = 13tt/40. 

spanning the interval < k < 400 £I p /va and 10 grid 
cells for the interval < < n/2. We assume cylindrical 
and reflectional symmetry. 

The wave energy per unit volume in k space is W k = 

{Bi-B k + E* k -[d(cog h )/dco]-E k }/(87z), where g h is 

the hermitian part of the dielectric tensor J^]. We 
evolve W k in time using "detailed energy conservation," 
i.e., by keeping track of the change in the particle en- 
ergy AS resulting from waves in each wavenumber 
bin, and deducting AS from the wave energy in that 
wavenumber bin. It can be shown that this method is 
equivalent to evolving the waves using the analytic for- 
mula for the damping or growth rate y k given by Hill for 
the limit in which \y k \ <C \(£> kr \. 

We integrate equation © for / (and the corre- 
sponding equation for W k resulting from detailed 
energy conservation) using an implicit time step- 
ping algortihm, the biconjugate gradient stabilized 
method lfl5ll . We hold the wave power spectrum W k 
fixed at the value 10" 7 fe"VBgv A /(8^) for 
= (13/40)?r. At all other values of 0, we initially 
set W k = l{r 10 Jfc -1 Vfi§v^/(8ffQ!) and then we allow 
W k to vary in time. The proton distribution function / 
is initially Maxwellian with a thermal speed v< v 2 > 
of 0.012 va- The minimum value of at the cell center 
in the run that we present is m j n = 7T/40. We note that 
the scattering contours for waves with = 0^,, are very 
similar to the 77 = constant lines. 

At early times, all the waves are damped by interact- 
ing with thermal particles as shown in Figure [3] Waves 
with = (13/40)?r, which are initially dominant, cause 
particles to diffuse along the relatively steep scattering 
contours shown by the dotted lines in Figure|4] When the 
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FIGURE 4. At t = 8,857a- 1 , the contours of constant / 
(solid lines) become steeper than the scattering contours of 
the waves with 6 = k/40 (dashed lines) and almost aligned 
with the scattering contours of the waves with 6 = (13/40)7T 
(dotted lines). Waves with 6 = k/40 become unstable and 
subsequently grow in amplitude. 



contours of constant / in the simulation become steeper 
than the scattering contours of the waves with 9 = m i n , 
waves with 9 = 9^ are amplified. After the energy in 
waves with 9 = m j n exceeds the energy in waves with 
9 = (13/40)7T, the small- 9 waves begin to dominate, and 
the contours of constant / start to align with the scatter- 
ing contours of the waves with 9 = m ; n , as shown in 
Figure [5] Even though we ongoingly input energy only 
into oblique waves, (quasi) parallel waves at 9 — 9 m i„ ul- 
timately dominate, and the distribution function evolves 
to a state in which / ~ /(fj) 

DISCUSSION 

Isenberg iflrjll has shown that A/IC waves with 9 = are 
unable to explain the heating and acceleration of protons 
in the fast solar wind, primarily because / relaxes to- 
wards a state in which / ~ /(t]), after which the pro- 
tons are only weakly heated by the waves. In contrast, 
we have shown that when protons are heated by oblique 
A/IC waves, the distribution function does not relax to- 
wards a state in which oblique-wave damping vanishes. 
Instead, / again approaches a state in which / = /(tj), 
and oblique waves continue to damp on the protons, 
causing protons to diffuse across 77 = constant surfaces 
in the v± — V11 plane. Because of this, oblique A/IC waves 
have the potential to be more effective than 0=0 waves 
at heating protons in the corona and solar wind. 
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FIGURE 5. At t = 14, 762^' , the amplitudes of the waves 
with 6 = k/40 exceed the fixed amplitudes of the waves at 6 = 
(13/40)tt and the contours of constant / (solid lines) become 
almost aligned with the scattering contours of the waves with 
9 = k/40 (dashed lines). 
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